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We propose a simple method for the deterministic genera- 
tion of an arbitrary continuous quantum state of the center-of- 
mass of an atom. The method's spatial resolution gradually 
increases with the interaction time with no apparent funda- 
mental limitations. Such de-Broglie Wave-Front Engineering 
of the atomic density can find applications in Atom Lithogra- 
phy, and we discuss possible implementations of our scheme 
in atomic beam experiments. 



Engineering of quantum states has been a widely dis- 
cussed topic for the last decade. Apart from a purely 
academic interest, there exist numerous applications of 
quantum state engineering including preparation of non- 
classical states of a cavity electromagnetic field, "pro- 
gramming" of a trapped- ion-based quantum computer, 
and atom lithography. Initial theoretical suggestions 
[|l]-[| for the preparation of a pre- chosen quantum state 
of a cavity field were based on the so-called conditional 
measurement method, where the target state is reached 
after a "successful" sequence of quantum measurements, 
while the "unsuccessful" measurement events are dis- 
carded. In the schemes |^J||l(| , applicable to both cavity 
light and external motion of a trapped ion, a two-level 
atom, coupled to the quantum field of interest as well 
as to a controllable external laser light, plays a role of 
a "bus", which transfers, in a prescribed way, popula- 
tion and coherence between the discrete eigenstates of 
the quantum field. Similar ideas were used to generate 
an arbitrary internal state of a multilevel atom Ac- 
cording to the suggestions |)],^| , the adiabatic population 
transfer process allows a one-to-one mapping between a 
quantum state of a Zeeman multiplet and a cavity field. 

The quantum state engineering methods listed above 
deal with systems of a discrete spectrum. In our pa- 
per we suggest a simple method to create an arbitrary 
continuous motional state of a free atom starting from 
a plane wave as an initial condition. The role of a bus, 
transferring the coherence between the initial and "tar- 
get" states, is played by an external uniform force field: 
the "target" motional state is encoded in the time depen- 
dence of the amplitude of applied laser light. 

The general idea of the set-up for the realization of 
our de-Broglie Wave-Front Engineering scheme was in- 
spired by the precision position measurement technique 
suggested and experimentally realized by J. E. Thomas 



[ pjj . Let us consider a two- level atom, interacting with a 
magnetic field H[z) = —az whose amplitude varies lin- 
early in space. Suppose that the internal atomic state 
|1) does not interact with the magnetic field (the corre- 
sponding Lande factor equals zero: g\ =0), whereas the 
Lande factor for the state |2) has a finite value gi — g. 
The energy difference between the states |2) and |1) will, 
thus, depend linearly on the position of the atom: 



^2,1(2) = ^2,1(0) - Tz , 



(i) 



where ^2,1 (0) = E2—E1 is the energy difference between 
1 2) and |1) in the absence of the magnetic field, and the 
gradient force acting on atom in state |2) is T = «/iBohr<7- 
In what follows we will assume, without loss of generality, 
that both a and g are positive numbers. 

Suppose for a moment that our goal is to create the 
narrowest possible position distribution of atoms in the 
state 1 2) centered at a position z = z*, and that the 
initial condition corresponds to a state |1) atom of some 
momentum po (whose value we can adjust at will). The 
simplest (but as we will see below not the optimal) way 
to approach the above goal is to apply a monochromatic 
spatially uniform laser field of a frequency uo = 6^2,1 (^*) 
for some period of time T (solid vertical arrow, inset for 
the Fig. |l|). As it is shown in the work jOJ, for a given 
value of the force T the spatial width of the distribution 
of atoms created in the state |2) is limited from below by 
a value 



d = 



2MT 



(2) 



(the so-called diffraction limit) no matter how long the 
interaction time T is. This limit is reached at a time of 
the order of 



[2%M\ 



(3) 



At times shorter than the time ([|), a wave packet 
of a minimum position-momentum uncertainty relation 
(Sz5p ~ Ti) is created. The external force T broadens 
the momentum distribution according to Sp ~ TT and, 
therefore, the spatial width of the state |2) distribution 
decreases with time as Sz ~ Tij TT . For long interaction 
times though, the peak width starts increasing quadrat- 
ically as the interaction time increases (see for example 
the spatial distribution calculated for the case T = 4.8r 
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shown at the Fig. [l]). Such a broadening is caused by both 
the quantum-mechanical diffraction of the wave-packet 
being prepared and the acceleration of the wave-packet. 
In what follows we will show that it is possible to suppress 
both diffraction and acceleration broadening using a sim- 
ple modification of the time dependence of the laser field 
amplitude (which was time-independent in the above ex- 
ample) . 

Let us write the field amplitude V(t) in a form 



V(t) = Ve 



-i J* uj2,i(z s {t' ,po,z*)) dt' 



(4) 



where the trajectory of the "resonant" point z s (t,po, z*) 
should be optimized in such a way that at the final time 
T the state |2) atoms will form a narrow peak, centered 
at z = z*. Notice that at a given time t the field (Q) 
plays a role of a spatially localized source of atoms in 
internal state |2) and with momentum Consider then 
a classical analog of our problem: 

Find a trajectory z s (t,po, z*) of a classical source 
of atoms of an initial momentum po 5 such that all 
the atoms emitted will reach the "target" z = z* 
at a preselected time T. Atoms are supposed to be 
affected by a force T. 

Such a trajectory does exist: it is given by 

Po(T-t) F(T-t) 2 



z a (t,po,z*) = z* 



(5) 



M 2M 

Let us now insert the ansatz (||) to the expression for 
the field amplitude (||) and evaluate the equations of mo- 
tion using this amplitude. At time T the state |2) com- 
ponent of the atomic wave function will be given by 



[i>2(z)]t=T = 




ip(z — z*) /h 



dpe 



e H( P0 +rT/2) z /K] gmc 



(6) 



[5z(T)\ ' 

where for all interaction times T, the spatial width of the 
distribution reads Sz(T) = 2H/FT, and no diffraction 
or acceleration broadening is present (solid line at Fig. 
p]). Here p- in is the initial density of atoms in the state 
|1). (Above we have omitted some insignificant constant 
phase factors depending on the definition of the internal 
states |1) and |2).) Note that the smc-shape is ultimately 
the best approximation for a ^-functional peak S(z — z*) 
one can create using a [po ; po + TT\ window in momen- 
tum space. 

The ansatz (^,^) motivates our strategy for "de-Broglie 
wave- front engineering" of motional quantum states. 
Imagine that one's goal is to prepare an atom in a mo- 
tional state <j){z) (normalized to unity: J dz \(j)\ 2 — 1). 
Let us represent the "target" state as a continuous super- 
position of the £-peaks: <p(z) = f^^dz* 5(z — dz*)(j)(z*). 
The state engineering process will involve then the fol- 
lowing steps: 



1. Prepare the atom in the internal state |1) and in an 
external state corresponding to the p = po eigen- 
state of the atomic momentum 



t=0 



/pi n e 



+ip z/h 



(7) 



where the initial momentum po is chosen in such 
a way that the momentum window [po! Po + TT\ 
covers entirely the momentum distribution 



Hp) 



i r +oc 



dze- ipz/h 4>(z) 



(8) 



of the "target" state, T and T being the typical 
magnetic field gradient and typical interaction time 
available in the given implementation; 

2. Apply for a time T a laser field 

V(t) = y^ e -iJ*«^Mt', Po ,o))dt' ^ (g) 

where 

V( t )= J^rL ^ dz* e -*< T *°>** x 



2tt^(0) 

^—i J*[uj2,i(z s (t,po,z*))—uj2,i(zs(t,po,0))] dt' 



VoFT 



x 4>(p(t,Po)) 



(10) 



is the field amplitude adjusted to the "target" state, 
and p(t,po) = Po + T(T - t); 

3. At the time T measure the internal state. If the 
atom remains in the state |1), repeat the above 
steps. If atom is detected in the state |2), the prepa- 
ration procedure is complete. 

After a lengthy but straightforward calculation one can 
show that in the course of the preparation procedure the 
initial state gets transformed to a state 



t=T 



f 



Po+FT 

dp (j){p) e +ipz/h 

PO 



(*), 



(11) 



very close to the "target" state 0, which was, we recall, 
assumed to be localized in momentum space within a 
range covered by the [po; Po + TT\ interval. This is the 
central result of our paper. 

Notice that the field amplitude @ is, apart from an 
overall time-independent amplitude, a product of two dis- 
tinct time-dependent factors. The second one is not spe- 
cific for a particular "target", but only for a given mo- 
mentum window [po; Po + TT\. This factor can be set 
once and for all for a broad class of targets. We re- 
call that the purpose of this factor is to compensate the 
quantum-mechanical diffraction and acceleration of the 
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distribution being generated. The first factor (10) cor- 
responds to a finite duration pulse whose amplitude is 
proportional to the target wave function in the mo- 
mentum representation: it allows encoding of the "tar- 
get" wave-function in the spectrum of the applied laser 
field. 
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FIG. 1. Spatial distribution of atoms in the state |2) 
(ab initio quantum- mechanical simulation). The desired pat- 
tern corresponds to a functional peak centered at z = 0. 
Dashed line is a result of an application of a monochro- 
matic field at resonance at z = 0. Solid line shows the 
"de-Broglie wave- front engineering" result, where the fre- 
quency was chirped according to "classical" ansatz (|^,^) taken 
at po = 0. Dashed-dotted line shows the initial distribution 
of the state |1) atoms. Interaction time is T = 4.8 x r, field 
amplitude is V = 1. x (h/r). The characteristic length d and 
time r of the problem are defined by the expressions (^|) and 
(^) respectively. The inset illustrates the basic idea of the 
"de-Broglie Wave-Front Engineering" method: the trajectory 
of the resonant point is designed in such a way that state |2) 
atoms created at different stages of the process come to the 
target at the same time t — T. 

Let us discuss now possible implementations of our 
scheme in the field of atom lithography. Consider a 
set-up depicted in Fig. |^a. We assume that a beam 
of atoms in state |1) with a mean longitudinal veloc- 
ity (v x ,o) interacts with a superposition of a Gaussian 
laser beam of a time-dependent amplitude and two mag- 
netic fields (z- and x-dependent respectively). The am- 
plitude of the x-dependent component of the magnetic 
field H2(x) = —(3(x — £det) 2 is chosen to be equal to 
P = A^Bohr go? /2M(^ 5 o) 2 • The z-dependent component 
of the magnetic field is the same as in the previous ex- 
ample. 

Consider first an idealized case of a monochromatic 
(v x ,o = ( v x,o)) ideally collimated (v Vy o = v z $ = 0) atomic 
beam, vo being the velocity of atoms in the incoming 
beam. Within the paraxial approximation and in the 



interaction picture with respect to the x-dependent Zee- 
man shift, the laser field "seen" by a moving atom reads 



V{t) = V{i + t i: 



(12) 



where hu 2 ,i(z) = 7kj 2 ,i(0) - Tz + R, R = h z k 2 /2M is 
the recoil energy, ki aser = ke y is the wave vector of the 
incident light, t = t — t[ n is the time counted from the 
"entering" time t- m when the atom enters the interaction 
zone [x- m < x < £det] 7 the left edge of the interaction zone 
is conventionally defined as x- in = x^et — 2(^det — ^laser), 
the interaction time is given by T = (xdet — X\n)/v x $, 
the "source trajectory" is defined by the expression (g), 
and the time width of the Gaussian pulse is Ti aser = 
/ {v x $) where w is the spatial waist of the Gaus- 
sian beam. 

Imagine that one's goal is to generate, by chosing a 
proper time dependence of the field amplitude V(t), a 
pre-chosen "target" density profile a(z) (scaled in such 
a way that cr(0) = 1) in the detector plane x = x^et 
(or, equivalently, at t = T). The quantum state engi- 
neering recipe pOj) does not directly apply: according to 
this recipe the function V(t) should be a finite-duration 
time pulse, and the probability for atom to join the in- 
teraction zone at the time of the pulse would vanish for 
long exposure times. In other words, the field amplitude 
( |l2|) contains a new (as compare to fl|l0|)) uncontrollable 
parameter t' m , randomly distributed within an interval 
[0; T], T being the "exposure time". We suggest now to 
replace the pulse jlC|)) by a periodic sequence of pulses. 

Let us consider the following ansatz for the field am- 
plitude: 



V(t) 



+N 

Vq ^2 V a ( z n)e~ 

n=-N 



(13) 



where z n = Azn, and oj n = c^2,i(^n); the space domain 
[— Az TV; +Az N] is supposed to cover the target's profile 
a{z). For long enough deposition times T ^> N/Alu the 
surface density of the state |2) deposited on the detection 
plane will be given by a sum of Gaussians profiles with 
an envelope given by the "target" pattern: 



n(y,z) « n 



(£-£nT 



% )e 2( 5 *)2 &n Q{5z/Az)o-(z) : 



-N 



where the width of an individual Gaussian is given by 
5z = y/Trh/FTiaser = ^/2h(v x j o) / ' Tw the overall den- 
sity amplitude is n = (jm^)(VTi aser /ft) 2 , j in is the 
flux density in the incoming atomic beam, and Q(() = 
^n'=-oo ex P[ — { n ') 2 I^C 2 ]' Notice that the produced pat- 
tern does not depend on particular values of the entering 
times ti n of the individual atoms, as a result of averag- 
ing over a uniform distribution of t- m within an interval 
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[0; T]. As in the previous example the spatial resolu- 
tion Sz of the atomic pattern engineering method is lim- 
ited only by the position-momentum uncertainty relation 

Sz ~ h/FTi^er. 

In Fig.|b we show the result of our attempt to generate 
the "boa swallowed an elephant" pattern , using this 
method. 

For a realistic atomic beam of a finite spread 
in the longitudinal (5v Xj q) and transverse (5v Zj q and 
Sv Vy o) velocities, the spatial resolution of the above 
lithographic method will be limited by S^ Vx0 ^z ~ 
(J r T 2 /M)(5v x ,o/(v x ,o}), fi(5v z , ) z ~ $ v z,oT, and 
s (Sv y ,o) z ~ k5vy y o/f, respectively. 
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FIG. 2. A possible implementation of the "de-Broglie 
Wave-Front Engineering" technique in atom lithography, (a) 
Schematics of the set-up. (b) Numerical simulation of an at- 
tempt to reproduce the "boa swallowed an elephant" pattern 
[12] using the de-Broglie engineering technique. Dashed line 
shows the spatial distribution of atoms in the state |2) in the 
detection plane. Interaction time is T = 1. X r, Rabi fre- 
quency is Vq = 1. X h/r. Time width of the Gaussian pulse is 
Tiaser = 0-3 x r. The pattern is represented by 2N + 1 = 127 
pixels on a Az = 3.7 x d grid. The distribution is a result 
of averaging over 20 Monte- Carlo realizations corresponding 
to the "entering time" £i n randomly distributed within an in- 
terval [0; 27r/Aa;], where Aou — TAz/fi. The characteristic 
length d and time r of the problem are defined by the expres- 
sions (0) and (pi) respectively. 



To conclude, we have presented a method for gener- 
ation of an arbitrary motional quantum state of a free 
atom. For a target state <j>(z) such that its momen- 
tum representation is well localized within an in- 
terval [pmin; Pmax], our de-Broglie wave- front engineering 
method allows one to create a state 

e i Po z/h _^ / dp4>{p) e ipz/h * (j){z) , 

where the initial momentum po and the interaction time 
T are chosen in such a way that the interval [p m in] Pmax] 
belongs to the [po; Po+3~T] window, T being the gradient 
force used in course of the generation procedure. 

We foresee that, using the classical- mechanical analogy 
similar to (||), our method can be easily generalized to the 
case of trapped particles. 

We present also a modification of the method, suitable 
for lithography with atomic beams. The spatial resolu- 
tion of the suggested lithographic method is limited only 
by the position-momentum uncertainty relation, and nei- 
ther quantum-mechanical diffraction of the pattern nor 
acceleration of atoms during the generation process af- 
fects the resolution. 

Let us give some realistic estimates for the spatial 
resolution one can achieve using the above lithographic 
method; we will use Fig. [l| as an example. Recall that 
this plot shows the narrowest peak, which could be ob- 
tained for given values of the field gradient T and in- 
teraction time T. For the Argon mass M = 30amu and 
a realistic value of the magnetic-field-induced (alterna- 
tively the Stark shift in a spatially varying laser field 
induced [jy]) gradient force T jU = 2tt x 10 9 Hz/cm (2tt x 
10 12 Hz/cm) the natural units of length and time will 
be given by d = (h 2 /2MJ r ) 1 ^ = 108.nm (10.8nm) and 
r = (2hM/J=" 2 ) 1/3 = 14.7/is(.147/is), respectively. The 
HWHM of the peak shown in Fig. |l| will correspond then 
to 63.0nm (6.3nm) obtained for an interaction time of 
T = 4.8r = 70.8/is(.708/is). Note also that lithography 
techniques sensitive to the internal state of atoms do ex- 
ist and they were experimentally demonstrated for the 
case of metastable atoms JT^f . 
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